Phase separation in supersolids 
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We study quantum phase transitions in the ground state of the two dimensional hard-core boson 
Hubbard Hamiltonian. Recent work on this and related models has suggested "supersolid" phases 
with simultaneous diagonal and off-diagonal long range order. We show numerically that, contrary 
to the generally held belief, the most commonly discussed "checkerboard" supersolid is thermo- 
dynamically unstable. Furthermore, this supersolid cannot be stabilized by next near neighbour 
interaction. We obtain the correct phase diagram using the Maxwell construction. We demonstrate 
the "striped" supersolid is thermodynamically stable and is separated from the superfluid phase by 
a continuous phase transition. 
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Nearly forty years ago, Penrose and Onsager [[[] posed 
the question: Is it possible for a bosonic system, like 
4 He, to have a phase where long range crystal order 
("solid") and off diagonal long range order (superfluid- 
ity) co-exist? Their answer, based on an analysis which 
ignored zero-point fluctuations, was that such supersolid 
phases do not occur. Subsequently, it was arg ued |-|] 
that including the effect of large zero-point quantum fluc- 
tuations in the crystal phase might allow for the exis- 
tence of a supersolid. This question has continued to 
spark much theoretical, numerical and experimental in- 
terest. H There has been a convergence of agreement, 
based on mean- field §,HH and numerical |§-[L2| work, 
that supersolids do exist, in 2c? lattice models, particu- 
larly in systems in which the density of bosons is doped 
away from the commensurate fillings which are optimal 
for charge ordering. The existence of supersolids for 2-d 
quantum bosons has fundamental implications to vortex 
phases in superconductors because of formal mappings 
between the problems. [ [l3| In this paper, we demonstrate 
that the most discussed of these lattice supersolid phases 
is thermodynamically unstable, and argue that it does 
not exist in any region of interaction strength or density. 

Consider a 2-d square lattice with one hard-core boson 
for every two sites (p — \) interacting with near neighbor 
(nn) repulsion. If the interactions are weak, the bosons 
will be mobile and condense into a superfluid phase at 
low temperature. If repulsion is strong, the system will 
freeze into a charge density wave pattern in which sites 
are alternately occupied and empty. At p = k these pos- 
sibilities are mutually exclusive. If we remove or add a 
boson, the resulting bosonic defect could "hop" among 
the background of charge ordered particles if the zero- 
point fluctuations are large enough. A dilute gas of such 
defects may Bose condense and form a superfluid super- 
imposed on the background of crystal order, a "super- 
solid" phase. If, instead, next near-neighbor (nnn) re- 
pulsion dominates, the charge ordering is in stripes, but 



the basic issue of a condensation of additional bosons 
coexisting with a striped pattern is as for checkerboard. 
While it is useful to think of separate frozen and super- 
fluid bosons, these quantum particles are indistinguish- 
able. All the bosons simultaneously participate in both 
types of long range order. 

Calculations supporting this intuitive physical picture 
are primarily based on mean field theory with spin wave 
stability analysis. They initially dealt with checkerboard 
charge order where the ordering vector for the structure 
factor is k = (tt, tt). Liu and Fisher Q argued that the 
supersolid phase exists for hard-core bosons with nn re- 
pulsion, but that it is unstable in the sense that the crit- 
ical velocity vanishes. If nnn repulsion is present, the 
supersolid can be stabilized. jl5| Numerical simulations 
of the quantum phase model (QPM), |l4| which describes 
soft-core bosons, showed that the (tt, tt) supersolid phase 
exists even without nnn repulsion, due to the soft cores. 
The supersolid is present even at half-filling, i.e. in the 
absence of any defects [fiof . Simulations of the hard-core 
bosonic Hubbard model similarly found that the super- 
solid phase exists in the absence of nnn repulsion off half- 
filling, but unlike the QPM is absent at half-filling. In 
addition, a mean field with spin wave analysis showed 
that this supersolid phase has a finite critical velocity. 
pl|,p[ However, it is generally accepted that at least in 
the presence of nnn repulsion, the (tt, tt) supersolid phase 
is stable. 

Discussions of stability based on nonvanishing critical 
velocity examine the effects of low energy excitations on 
an existing supersolid phase. To our knowledge, how- 
ever, there has been no discussion or numerical verifi- 
cation of the underlying thermodynamic stability of ei- 
ther of these supersolid phases against phase separation. 
Simulations done at fixed particle number which found si- 
multaneous diagonal and off-diagonal long range order, 
||[l0 11 1^,0, do not address the possibility of phase 
separation. In what follows we examine thermodynamic 
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stability of the checkerboard and striped supersolids by 
constructing the chemical potential-particle number re- 
lation and calculating the compressibility. 

We use a new Dual Quantum Monte Carlo Algorithm 
[pll to simulate the hard-core bosonic Hubbard model, 
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di (aj) are destruction (creation) operators of hard-core 
bosons on site i of a 2-d square lattice, and rij is the 
density at site i. The hopping parameter is chosen to be 
t = 1 to fix the energy scale. V\ (V2) is the near neigh- 
bor (next near neighbor) interaction. At V2 = 0, and 
after an appropriate sublattice spin rotation, this boson 
model is equivalent to the spin-^ antiferromagnetic XXZ 
model. In this language, superfiuid order corresponds to 
magnetic order in the XY plane, while density order cor- 
responds to magnetic order in the Z direction. 
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FIG. 1. The structure factor, .S^vr, 7r), and p s as a function 
of fixed density. The half-filled point p = 0.5 is a solid with 
p s — 0. For p close to p = 0.5, S(tt, tt) and p a are both 
nonzero. Moving away from half-filling, eventually S(tt, tt) 
will no longer scale linearly with system size and the system 
is a pure superfiuid. 

Traditionally, to determine numerically the nature of 
the ground state of (1), we evaluate, at fixed density, the 
superfiuid density, p s , and the equal time structure factor 
at the ordering vector q, 



(2) 



Ground state results for S(q) and p s are shown in Fig. 1 
(Fig. 2) for the checkerboard (striped) phase. In both 



cases p s is nonzero everywhere except precisely at half- 
filling, but S(q) also remains large off half-filling, indi- 
cating solid order. Using finite size scaling to extrapolate 
to the limit of infinite lattice size for these fixed density 
systems, ||,[l^,|ll],[il|[l7| one can show that density corre- 
lations arc indeed still long ranged off half-filling where 
p s 0. The conclusion is that checkerboard and striped 
supersolid phases exist in the thermodynamic limit. 

Already, however, it was remarked [|| that the energy 
versus density curves had small negative curvature in 
the (tt, tt) supersolid phase. It was speculated that this 
was evidence for phase separation. The (n, 0) supersolid 
phase showed no such negative curvature. It was recently 
shown numerically jl7) that the easy-axis spin- 1/2 XXZ 
model on a square lattice exhibits a first order spin-flop 
transition. These results confirm the discontinuous na- 
ture of the transition from the superfiuid phase as p is 
adjusted in the absence of nnn repulsion. 
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FIG. 2. The structure factor S(0, tt) and p s as a function 
of fixed density. The half-filled point p — 0.5 is a solid with 
p s — 0. For p close to p — 0.5, S(0, tt) and p 3 are both 
nonzero. As with the checkerboard case, finite size scaling for 
S(0, tt) determines the density at which solid order vanishes. 

To address the possibility of phase separation sys- 
tematically, i.e. to obtain the phase diagram in the 
interaction-chemical potential (p) plane, we must ob- 
tain p as a function of p, and use the Maxwell construc- 
tion ]1S[| . We first study the checkerboard case by fixing 
V% = and scanning the filling for several values of V\. 
The chemical potential for n bosons is calculated from 
the total energy: p(n) = E(n + 1) — E(n). We work 
on lattices up to size 12x12, and temperatures as low as 
j3 = 6 to access the ground state properties. 

Fig. 3 shows p versus p for V\ = 3, V% = 0, as in Fig. 1. 
The slope of this curve is the compressibility, n = dp/ dp. 
Two k < branches are clearly seen just before and after 
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the energy gap. The gap itself corresponds to the incom- 
pressible (w, 7r) solid at half filling and seen in Fig. 1. 
Using the Maxwell construction we find the critical value 
of the chemical potential, p c (vertical dashed line), and 
read off the critical filling p c in Fig. 3. The structure 
factor (Fig. 1) begins a very rapid rise at the point where 
k turns negative. ]2(| It is crucial to note that p s and 
S(tt,w) are both non-zero, Fig. 1, only in the unstable 
k < region (Fig. 3). On one side of the transition the 
phase is purely superfluid while on the other side it is 
a (tt, tt) solid. What was previously accepted to be the 
checkerboard supersolid at fixed density, lies entirely on 
the k < branches and therefore phase separates into a 
mixture of solid and superfluid phases at densities p = p c 
and p = \- The metastable states in Fig. 3 correspond 
either to the superfluid or the gapped insulating phases. 
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FIG. 3. p versus p showing the k < regions. The vertical 
dashed line shows the location of the transitions from the ^ 
Maxwell construction. Inset: The phase diagram for V% = 0. 
The solid line shows the continuous transition to the Mott 
phase at full filling, the dashed line shows the discontinuous 
transitions from the superfluid to the checkerboard solid at 
half filling. The density changes discontinuously across this 
line. The tip of the lobe is a continuous critical point. 

To check if nnn repulsion stabilizes this phase against 
phase separation, we did simulations [p(if with Vi = 3 
and V2 ranging from deep in the (tt, it) solid region to 
close to the boundary with striped order. We found the 
same k < behavior as in Fig. 3. Next near neighbour 
repulsion does not stabilize the checkerboard supersolid 
phase. 

Repeating the simulations that gave Fig. 3 for differ- 
ent values of V\ we construct the phase diagram in the 
(p/Vi,t/Vi, V2 = 0) plane, shown as the inset to Fig. 3. 
As the tip of the lobe is approached, the energy gap opens 
without k < regions in the p, p plane. Therefore this 



point is apparently a continuous transition. This is con- 
sistent with Ref. while Ref. [^2| finds a first order 
transition in a model with longer range (Coulomb) inter- 
actions. 

The same analysis for V\ = 0, scanning V% and p, de- 
termines the stability of the striped supersolid. Fig. 4 
is a typical plot of p versus p traversing the incompress- 
ible (gapped) striped solid at p = \. This is strikingly 
different from Fig. 3. There is no k < region: The 
phase transitions are all continuous. Furthermore, as p 
is increased from the lowest shown value, the slope, i.e. 
k, changes markedly at (p w 0.74, p w 0.25). We find 
that S(0, 7r) (Fig. 2) begins a rapid increase, indicating 
long range striped order, at precisely this particle den- 
sity, p. Since the superfluid density, p Sl is still finite, we 
conclude that the rapid crossover in k, like the behavior 
of the structure factor, signals the continuous transition 
from the superfluid to the (0,7r) supersolid. Increasing p 
further takes the system into the gapped (0, 7r) insulator, 
then back to the striped supersolid and finally to the su- 
perfluid phase at p « 19.2. At strong coupling, the gap 
(the jump in p) across the incompressible striped phase is 
4V2 = 20 for the 2-d square lattice. This value is reduced 
by quantum fluctuations as V2 decreases, eventually dis- 
appearing entirely at weak coupling. The absence of neg- 
ative compressibility regions indicates that all phases, in 
particular the striped supersolid, are thermodynamically 
stable. 
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FIG. 4. p versus p. Inset: The phase diagram for Vi = 0. 
The narrow regions sandwiched between SF and (n, 0) solid 
phases are the stable supersolid phases. 

Repeating these simulations for various values of V% 
gives us the phase diagram in the (p/V2, t/V2, V\ = 0) 
plane (inset to Fig. 4). The regions sandwiched between 
the superfluid and the striped solid are the two stable 
striped supersolid phases. As the tip is approached, the 
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supersolid phase gets narrower since the SS-(0, 7r)-solid 
transition approaches the SS-SF transition. This pre- 
vents us from resolving these transition points near the 
tip. The difficulty in the numerical determination of 
these points stems from the fact that when the SF-SS 
and SS-striped solid transitions get very close to each 
other, they start behaving numerically as multicritical 
points. However, it appears that the supersolid phase 
completely surrounds the striped solid phase except at 
the tip and the base where we have multicritical points. 

We also studied the effect of nn repulsion on the (0, tt) 
phase and found that the (0, n) supersolid remains sta- 
ble and that additional gapped phases appear at other 
special fillings. 1 10 2l| Whether there are associated su- 
persolid phases is under investigation [^3| . 

The fact that it is easier to support nonzero p s in a 
striped solid than a checkerboard one can be qualita- 
tively argued as follows: In a striped solid doped away 
from half-filling, defects have channels in which they can 
move at no interaction energy cost, and, importantly, the 
kinetic energy of these defects is set by t and can be con- 
trolled independently of the strength of the interaction 
V-i which determines the solid order. In a checkerboard 
solid, the motion of a defect proceeds through an inter- 
mediate state of energy 2Vi, giving a reduced effective 
hopping t G ff w i 2 /2Vi. V\ controls simultaneously the 
defect kinetic energy and the tendency to charge order. 
As a consequence, there is reduced ability to tune to a 
supersolid phase. It is still remarkable, though, that the 
striped phase forms even at very low densities. Indeed, 
in the fermion Hubbard model, very small doping (just 
a few percent) away from half-filling destroys long range 
spin order (antiferromagnetism), leaving little possibility 
that it might coexist with superconductivity. 

In this paper, we have presented Quantum Monte 
Carlo results for ground state correlations in the hard- 
core bosonic Hubbard model with near- and next near- 
neighbor repulsion. We show that the q = (tt, tt) checker- 
board supersolid, contrary to current beliefs, is an un- 
stable phase and does not exist thermodynamically for 
this model for any filling or nnn repulsion. Instead, the 
system phase separates into solid and superfluid phases. 
This contradicts mean-field predictions which examine 
stability via an evaluation of the critical velocity for 
spin waves, and find that in the presence of next near 
neighbour repulsion the supersolid is stable. We have 
not examined the soft-core case in detail, but prelimi- 
nary results indicate negative compressibility regions in 
that case too. [|| The quantum phase model for soft core 
bosons also exhibits negative compressibility phases. |24|] 

We found the striped supersolid phase, q = (tt, 0) to be 
stable and separated from the superfluid phase by a con- 
tinuous transition. The energy E(n) provides a signal 
of the transition: The compressibility exhibits a rapid 
cross-over from the superfluid to the supersolid phase, 
with ksf < kss- The issue of the stability of possi- 



ble supersolid phases at other densities and wavevectors 
which are associated with the presence of long range in- 
teractions is a fascinating one which is presently under 
investigation. 
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